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Abstract 

The nuclear quadrupole resonance (NQR) spectrum of strontium doped 
La2Cu04 surprisingly resembles the NQR spectrum of La2Cu04 doped with 
excess oxygen, both spectra being dominated by a main peak and one principal 
satellite peak at similar frequencies. Using first-principles cluster calculations 
this is investigated here by calculating the electric field gradient (EFG) at the 
central copper site of the cluster after replacing a lanthanum atom in the clus- 
ter with a strontium atom or adding an interstitial oxygen to the cluster. In 
each case the EFG was increased by approximately 10 % leading unexpectedly 
to the explanation that the NQR spectra are only accidentally similar and the 
origins are quite different. Additionally the widths of the peaks in the NQR 
spectra are explained by the different EFG of copper centres remote from the 
impurity. A model, based on holes moving rapidly across the planar oxygen 
atoms, is proposed to explain the observed increase in frequency of both the 
main and satellite peaks in the NQR spectrum as the doping concentration is 
increased. 
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A large amount of information about static and dynamic properties of high-temperature 
superconducting oxides has been provided by nuclear magnetic resonance (NMR) and nu- 
clear quadrupole resonance (NQR) experimentsB. In certain materials, however, some NMR 
and NQR results are still not completely understood. In part, this is related to the lack of 
theoretical support since, to date, only a few first-principles calculations have addressed the 
local charge and spin density distributions. An example of these unresolved problems is the 
interpretation of the NQR spectra in doped La2Cu04. This material becomes superconduct- 
ing only upon alloying on the La site or changing the oxygen content. There is therefore 
considerable interest in understanding the changes in the local electronic structure that the 
various dopants induce. 

Experimental observations! of NQR spectra for stoichiometric La2Cu04 show one sharp 
peak at 33 MHz assigned to ^'^Cu. Upon doping with strontium, a second weaker hue appears 
at a somewhat higher frequencylHi. This was called B-line to distinguish it from the main 
resonance line (A). Increasing the concentration of dopants causes both peaks to broaden 
as well as increase in frequency slightly. Similar behaviour is observedili when oxygen is 
used as the dopant. This produces an additional line which is so close in frequency to the 
B-line frequency of the Sr-doped spectra that both lines were attributed to the same physical 
origin. Hence it was proposed that holes are localized in CuOe octahedra adjacent to the 
charged dopants. These holes, if pinned at least for times of the order of microseconds, 
would then cause the NQR B-lines. This interpretation was corroborated by results from 
ab-initio calculations^ on the magnitude of the copper NQR frequency which showed that 
the Sr- and 0-dopants themselves had very different effects on the NQR frequency but 
the frequency shift for a copper adjoining a site occupied by a localized hole matched the 
experimental results very closely. Although there is strong evidence for various kinds of 
charge inhomogeneities in doped lanthanum cuprate, the occurrence of static pinned holes 
proved elusive. In addition, there are some subtle differences between the NQR spectra of 
Sr- and 0-doped samples, such as the asymmetry of the electric field gradient of the B-lines, 
which are not explained by the above model. 



In this work, we report on results of density functional cluster calculations that address 
the changes in the electric field gradients (EFGs) in La2Cu04 upon doping. Four different 
cluster models were used to simulate various experimental situations: 

(Ma) Cu06/Cu4Laio represents the La2Cu04 crystal and used to obtain the EFG at the 
central copper site 

(M/3) Cu06/Cu4La9Sr represents the La2-xSra;Cu04 crystal in the dilute limit x 
(M^) Cu06/Cu4Laio and Cu06/Cu4La9Sr, as above, but where positive charges are placed 
on each planar oxygen. This model is proposed to investigate the influence of holes that are 
moving rapidly across the planar oxygen sites 

{Ms) Cu07/Cu4Laio, the same as Cu06/Cu4Laio but where an extra oxygen atom occupies 
an interstitial site, representing the La2Cu04+5 crystal in the dilute limit 5 — > 0. 

Recently, extensive cluster studies for the pure La2Cu04 system have been presente d0. 
Spin-polarized calculations at the Hartree-Fock (HF) level and with the density functional 
(DF) method with gradient corrections to the exchange and correlation functionals have 
been performed for clusters comprising up to 9 Cu atoms in a Cu02-plane. It has been 
shown that an acceptable estimate of the EFG at the copper site can already be obtained 
with DF calculations with a small cluster (Cu06/Cu4Laio) where a (Cu06)^°^-ion is treated 
with an all-electron basis set and 4 surrounding Cu and 10 La atoms are represented by 
bare pseudopotentials. This is the cluster which is chosen as model in the present work. 
It is embedded in a lattice of 2464 point charges. The positions of some point charges were 
changed in such a way that the correct Madelung potential in the central region of the 
cluster was reproduced. All atomic positions were located according to the tetragonal phase 
with lattice constants0 a = 3.77 A and c = 13.18 A and with a Cu-O(apex) distance of 
2.40 A. For the cluster atoms standard triple-zeta (6-311 G) basis sets were employed. All 
calculations were spin-polarized and were performed with the Gaussian94 and Gaussian98 
programsEl. In the DF calculations the correlational functional of Lee, Yang, and ParJil 
was combined with the Beckelil exchange functional, which also includes corrections based 
on the gradient of the electron density. 



For model the density functional calculations with the BLYP functional give for Vzz-, 
the z-component of the EFG, at the Cu site a value of 1.419 a.u.. This corresponds to a ^'^Cu 
quadrupole frequency of 35 MHz if a quadrupole moment Qi^^Cn) = —.211 b is adopted^. 
The experimental ^^Cu NMR spectrum of pure La2Cu04 showsi a narrow peak at 33 MHz. 
The good agreement between theory and experiment may be accidental in view of the small 
cluster size. Furthermore, there is no general consensus about the precise value of Qi^^Cn). 
Other sources of inaccuracy are due to the limited basis set used for the representation of the 
atomic orbitals and the choice of the exchange and correlation functionals for the treatment 
of the generalized gradient approximations beyond the local density approximations. Since 
we focus in this paper on the relative changes of the local electronic structure around the 
copper atom which are induced by models for various doping effects, all calculations are 
carried out at the same level of approximation. We expect then that the difference between 
EFG values calculated for different clusters will not significantly be influenced by the above 
inaccuracies. 

It should further be remarked that the EFG cannot be explained in a simple one-electron 
(or rather one-hole) picture. The theoretical value results from cancellations between rela- 
tively large individual contributions from the various atomic shells which are listed in the 
second column of Table |l[ 

It has been shown0 that other choices of the exchange-correlation functional have little 
influence on Vzz- HF calculations, however, produce a T^^- value that is 36 % higher. This 
discrepancy between DF and HF results for the EFG at Cu sites has also been found in 
previous calculations for the YBa2Cu307 system. It is further discussed in Ref.llUli but is 
not the topic of the present work which focusses on the relative changes in Vzz upon doping. 

To simulate the influence of doping with Sr, the calculations were then repeated for the 
(Cu06)^°~-ion with the La'^"^ on top of the Cu replaced by Sr^^ (model M^). This simulates 
a situation where the hole introduced has moved to a place far from the cluster region or is 
distributed over several sites. The resulting changes in the contributions to Vzz with respect 
to the undoped system are given in the third column of Table |. The EFG is increased by 
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10.3 % which is very close to the observed difference of about 3 MHz in ^^uq between the 
additional peak (B-line) and the main peak (A-line). It should be noted that in models Mq, 
and M/3 the same atomic positions were used. Structural changes induced by replacing the 
lanthanum atom by strontium are thus not accounted for. 

At the Cu sites adjacent to the central Cu ion of the cluster, the replacement of La^^ by 
Sr^"*" leads to an increase of Vzz by 1 %. This explains the observed broadening of the main 
peak (A-line) in La2_aSr^Cu04. 

Similar cluster models (M^ for pure and Mp for Sr-doped La2Cu04) have previously been 
investigated by Martin0 with HF calculations with results which drastically differ from those 
of our DF calculations. We therefore also performed HF calculations for the clusters M„ 
and Mg. For the undoped system, our HF result {V^z = 1-94 a.u.) is in broad agreementlli 
with his {Vzz = 1-84 a.u.) and also with the value reported by Sulaiman et al.0 {Vzz = 2.23 
a.u.). For the Sr-doped system, however, the results differ strongly. In Ref.0, an increase of 
Vzz by 1.5 % was reported while our HF calculations predict a change by 5.9 %. 

Our DF results thus suggest that the B-line observed in La2Cu04 upon doping with 
strontium simply originates from the Cu atoms next to the positions of the dopant atoms. 
If the models Ma and Afg adequately describe the local electronic structure around the 
strontium dopant, a separate oxygen NQR spectrum is to be expected due to ^''O nuclei that 
bridge the strontium atom and its nearest copper atom. For model M^, our calculations 
give Vzz = 0.28 a.u. for the apical oxygens. For model they predict Vzz = 0.15 a.u. 
for the bridging 0(a). The number of oxygen nuclei that exhibit this spectrum, however, 
only amounts to pox/4: per formula unit where po is the abundance of ^''O nuclei in the 
sample. It should be noted, however, that the clusters used in this work were not designed 
for studying the EFGs at the oxygen sites. Indeed, using larger clusters for the undoped 
system, the values of the EFGs at the apical oxygens turn out to be smaller and closer to the 
experimental values of 0.22 a.u.. The pronounced calculated reduction of Vzz upon replacing 
a La by a Sr atom, however, is related to the charge redistribution and will persist also in 
more elaborate calculations. 
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Next we investigate the changes induced by an interstitial oxygen atom 0(i) at the 
location (a/2,0,c/4). In the tetragonal structure, this site is surrounded by four nearest 
oxygen neighbours at a distance of 2.09 A for the geometry parameters used here, followed 
by four La and four Cu. The latter form two pairs which lie in adjacent Cu02-planes and 
have the four 0(a) atoms at corresponding apical positions. This interstitial location for 
0(i) is undisputedil'ii. It is also evident that the squeezing in of the excess oxygen into 
the interstice will cause a relaxation of the atoms in the immediate neighborhood. The La 
positions seem to be unchangedil. The EFG value at the copper site depends crucially on 
the precise position of the apical oxygen near the interstitial site, even though the extent of 
its displacement is disputec&B. Therefore we determined this position of the 0(a) which 
also minimised the energy of our cluster model M^. For this, the cluster Cu06/Cu4Laio 
was extended by including an additional oxygen 0(i). This model also accounts for that 
situation where no additional holes are present in the cluster region. 

The positions of all cluster atoms and point charges of the background were first frozen 
except those of the apical oxygen with coordinates (x^, Ua-, Za) nearest to the 0(i) which was 
allowed to move on the plane Ua = 0- An energy minimum was found for Xa = —0.64 A 
and Za = 2.26 A. At this position, the distance to the interstitial oxygen is 2.73 A. This 
corresponds to a tilt of the 0(a) by 16*^ which is reasonable when compared to experimental 
information&B. Next the positions of the three other apical oxygens coordinated with 0(i) 
(which are represented by point charges only) were displaced accordingly. With this geom- 
etry, the EFG is no longer axially symmetric and the eigenvector of the largest component 
is rotated by 4° from the z-axis. The asymmetry of the EFG tensor amounts to rj = 0.05 
in approximate agreement with the experimental! value of 0.07. The largest eigenvalue of 
the EFG is 1.575 a.u. which is 11.0 % larger than in model M^- This corresponds to a 
shift in frequency of 3.8 MHz. HF calculations by MartinS which were carried out without 
taking the local distortions around the interstitial oxygens into account, showed an increase 
of about 5 MHz with respect to the results for the undoped case. 

Our model Ms can thus account for the occurrence and position of the satellite peak 

6 



in the NQR spectrum of La2Cu04+5. The close agreement of the frequencies of the B- 
hnes in strontium- and oxygen-doped lanthanum cuprate is therefore only accidental. The 
physical origins of the B-lines are quite different in the two cases. This is also reflected in the 
asymmetry parameter rj which is zero in model M/3 by symmetry. Repeating the calculations 
for a cluster with atomic positions chosen according to the orthorhombic structure yields an 
asymmetry of 77 = 0.01 for model Afg. 

Since each interstitial oxygen has four nearest copper atoms, the intensity of the B-line 
in La2Cu04+5 is thus expected to increase with 46 which is also in agreement with the 
experimental. 

Finally we study the shift of the A-line upon doping. We adopted model where holes 
rapidly move across the planar oxygen sites. For a rough estimation, positive point charges 
x/2 were placed near the planar oxygen sites in the cluster Cu06/Cu4Laio. To be precise, 
point charges with values x/8 were put at positions where the 2p^ and 2py orbitals of the 
planar oxygens attain their extremal values. The calculations were performed for x = 0.01, 
0.02, and 0.04, respectively. In the fourth column of Table |, the corresponding changes in 
the contributions to V^z with respect to the undoped system, which turned out to be linear 
in X, are listed for x = 0.04. We obtain a shift d{logVzz)/dx ~ 0.7 which is in reasonable 
agreement with the experimental value d{\og uq) / dx ~ 0.6 since no account is given to the 
fact that the lattice constants also slightly change with the doping concentration. Applying 
the same procedure to model Afg yields d(\ogVzz)/dx ~ 0.6. 

Our conclusions are as follows. The results obtained with the cluster Cu06/Cu4La9Sr 
(model Afg) predict an increase of the EFG at the Cu-site immediately adjacent to the Sr by 
10 %. This would explain the occurrence and position of the NQR B-line in La2-xSra.Cu04 
for X ^ 0. In addition, the observed broadening of the A-line with increasing doping 
concentration could be attributed to the small shifts calculated for Cu sites not immediately 
adjacent to the Sr. Next, the shift of the A-line with increasing x could be explained by holes 
moving rapidly across the planar oxygens (model M^). We have shown that by introducing 
excess oxygen in La2Cu044-5 (model Ms), the EFG at the Cu is also increased by an amount 
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which accidentally is about the same as doping with Sr produces. Hence our results disagree 
with previous interpretations of the Cu NQR spectra in the doped La2Cu04 system, in 
particular that of Martin0 which attribute the origin of the B-lines to localized holes that 
are introduced either by Sr or O doping. 

We express our gratitude to T. A. Claxton for enlightening discussions and critical read- 
ing of the manuscript. We gratefully acknowledge the help of P. Hiisser and E. P. Stoll. This 
work is partially supported by the Swiss National Science Foundation. 
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TABLES 

TABLE I. Contributions to the EFG component Vzz in atomic units from the nuclei within the 

cluster, from the point charges around the cluster, and from the individual shells. The "remainder" 
lists just the (small) rest of all other contributions that cannot be assigned to a particular shell. 
The second column shows the values for pure La2Cu04 (model Mq). In the third and forth column 
the changes are given which result when one La^'^ is replaced by Sr^+ (model M^) and for mobile 
holes corresponding to x = 0.04 (model My), respectively. 





pure 


Sr doped 


holes x=0.04 


Nuclei 


0.375 


0.003 


0.000 


Point charges 


0.012 


0.000 


-0.002 


Px) Py 


-2.393 


0.012 


0.013 


Pz 


0.770 


0.028 


0.003 


•^X^— J/2, dxy 


-15.163 


0.043 


0.011 




8.829 


0.040 


0.011 


dj/z) (izx 


9.003 


0.002 


0.001 


Remainder 


-0.012 


0.019 


0.007 


Total 


1.419 


0.147 


0.044 
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